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Abstract A method is presented that allows to solve the Faddeev integral equations of the 
semirelativistic constituent quark model. In such a model the quark-quark interaction is 
modeled by a infinitely rising confining potential and the kinetic energy is taken in a rela- 
tivistic form. We solve the integral equations in Coulomb-Sturmian basis. This basis facili- 
tate an exact treatment of the confining potentials. 
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1 Introduction 

The Faddeev approach to quantum-mechanical three-body problems has been efficiently 
applied to a variety of problems for many years. The main advantage over other methods 
consists of solving general Hamiltonians with two and three-body forces, implementing a 
priori the proper boundary conditions (bound state and/or scattering), and in observing the 
proper symmetries for (non)identical particles. While problems with short-range forces can 
be treated in a straightforward manner, special care must be exerted if long-range interac- 
tions are present. This applies in particular to the three-body Coulomb problem, but also to 
the case when confining interactions are present. The reason is that any long-range force 
influences the asymptotic behavior such that the system never becomes free of interactions. 
As a result, in the Coulomb case one has infinitely many bound states accumulating at zero 
energy, whereas for the confinement there exist infinitely many bound states with arbitrarily 
immense binding energies. 

Faddeev equations have already been applied in the past to confined three-quark sys- 
tems. The approach was first developed for a confinement of harmonic oscillator type, and 
one has achieved an exact solution for the three-body ground state 1 1 ]. The same method was 
subsequently followed for a non-relativistic three-quark Hamiltonian, including the addition 
of a linear confinement to a hyperfine interaction [213J. The technique always consisted of 
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defining the zero point of the spectrum such that the three-body states in question become 
deeply bound and standard methods become applicable. One of the notable disadvantages 
of this approach consists of difficulties in the numerical treatment as the partial-wave expan- 
sion becomes only slowly convergent. 

Some time ago, one of us proposed an alternative method for dealing with confine- 
ment in three-quark systems |4]. It is similar to the one followed in the Coulomb case 
and consists of splitting the quark-quark potential into a long-range confining part and a 
short-range (non-confining) interaction. The Faddeev procedure is then applied only to the 
non-confining short-range parts. The resulting integral equations were solved by using the 
Coulomb-Sturmian separable expansion approach. Subsequently the method was success- 
fully applied to a semirelativistic constituent-quark model |5|. 

In this approach one requires the evaluation of the matrix elements in the Green's oper- 
ator corresponding to the long-range part of the Hamiltonian between the basis states of the 
separable expansion. Whereas in the Coulomb case it is readily possible to compute these 
matrix elements in the Coulomb-Sturmian basis. An analogous result has not been known 
for the confining interactions, rather the confinement has been treated in an approximate 
way [4.5 1. The comparison with results for baryon masses obtained from variational calcu- 
lations suggested that this might not be a bad approximation, but one certainly desires to 
get rid of it, especially with regard to getting reliable wave functions for baryon ground and 
excited states. 

In this paper we present a method to solve the Faddeev equations for confining interac- 
tions without resorting to any such approximation. The technique makes use of a recent re- 
sult obtained for the Green's operator of a Hamiltonian consisting of interactions expressed 
as polynomials of the inter-particle separations (6J. In the following Section we shortly 
outline the Faddeev approach to three-body problems. Then, in Section [3] we explain the 
procedure for modifying the Faddeev approach in order to treat confining interactions. Sub- 
sequently we present the solution of the Faddeev equations through the Coulomb-Sturmian 
expansion method. We give the results for the nucleon and A ground and excited states, in 
the case of the relativistic Goldstone-boson-exchange constituent quark model |7). 



2 Faddeev approach to three-particle problems 

We consider a three-particle Hamiltonian 



where is the kinetic-energy operator and v a , cc = 1,2,3, are the mutual interactions of 
the particles. We represent it through the usual configuration-space Jacobi coordinates: e.g. 
xi is the coordinate between particles 2 and 3 and y t is the coordinate between the center of 
mass of the pair (2,3) and particle 1. Let us assume for a moment that the interactions v a 
are short-range potentials. 

We want to solve the Schrodinger equation 



H = H i0) + Vl +v 2 + v 3 , 



(1) 



#|f) = (# (0) +vi+V2 + v 3 )| l f')=£| l f') . 



(2) 



Introducing 



G (0) (£) = (Z 7_ // (0) ) -l 



(3) 
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and rearranging pi, we have 

|f> =G ( °)(£)(v 1+ v 2 + V3)|f> 

= GM(E)vi\V) +gW(E)v 2 \*F ) +G(°'(£)v 3 | , f>. (4) 

This means that the three-particle wave function f) naturally breaks down to three compo- 
nents 

1^) = l^l) + l^2> + |^3) , (5) 

where 

\y a ) = G®(E)v a \Y), a = 1,2,3, (6) 

are the Faddeev components. The Faddeev components satisfy the set of equations, the Fad- 
deev equations, 

{E-H® -vi)\Vi) = vi(|v^2> + |v^» (7) 
{E-H® -v 2 )\W2) =v 2 (|Vt) + |V3» (8) 

(£_ ff (Q)_ V3 )| ¥ 3) =V3(|Vl) + |V^» • (9) 

By adding up these equations and considering Q we recover the Schrodinger equation. 
So, the Faddeev procedure is no more and no less than a way of solving the Schrodinger 
equation. With the help of channel Green's operators 

G a (E) = (E-H^-v a )-\ (10) 



we can also recast the Faddeev equations into an integral equation form 



Vi)=Gi(£)vi(|YO + |Y*» (11) 
Wz) = G 2 (£)v2(IVi) + IV3» (12) 
V/ 3 > = G3(E)v 3 (|vi} + |V2» • (13) 



In general, the potentials may support bound states. So, we may have a situation when 
particles 2 and 3 form a bound state and particle 1 is far away. Alternatively particles 1 
and 2 form a bound state and particle 3 is away, and so on. So, in the three-particle wave 
function we may have three genuinely different two-body asymptotic channels. Although 
these two-body channels represent genuinely different physical situations, yet in solving the 
Schrodinger equation, we need to impose all these different boundary conditions on a single 
wave function. 

The Faddeev components posses a simpler structure. For example, in Eq. j6|, the short 
range potential vi acting on |f), suppresses those asymptotic structures when particles 2 and 
3 are far away, i.e. it eliminates the two-body channels when either particle 2 or particle 3 is 
at infinity. Consequently, | y/i ) contains only one kind of asymptotic channels, channels when 
particle 1 is at infinity. A similar statement is valid for i// 2 ) and \jfj). So, with the Faddeev 
decomposition we achieve an asymptotic filtering, i.e. we split the wave function into parts 
such that each component possess only one kind of asymptotic behavior. Then solving the 
Faddeev differential equation we need to impose simpler conditions on the components. 
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3 Faddeev equations for confining potentials 

It is obvious that if we apply the Faddeev decomposition to long-range confining potentials, 
we run into trouble. If v a in Eq. |6| is of long-range type, then applying on f), it would 
not suppress the components where particles /3 or y are far away. The procedure does not 
lead to any simplification in the asymptotic behavior of the components: the whole Faddeev 
procedure does not make much sense. 

A similar situation occurred when the Coulomb potential was directly inserted into the 
Faddeev equations. This approach turned out to be a complete failure. The way out of the 
trouble was proposed by Merkuriev JS). He proposed to split the Coulomb potential into 
long-range and short-range parts and applied the Faddeev procedure only to the short range 
parts. The set of Faddeev-Merkuriev integral equations possesses a compact kernel and its 
asymptotic analysis provided the boundary conditions to the differential equations |8|. The 
three-particle systems with attractive Coulomb potentials became amenable to numerical 
treatment both in differential |9, 10 1 and integral equation formalism ifTTl . 

Inspired by Merkuriev's treatment of the Coulomb potential we adopted [4| a separation 
of the confining potential into long range confining and short range non-confining terms 

Va = Va ] +V ( a\ (14) 

where v a ' and v a ' are confining and short range potentials, respectively. Then we can write 
the Hamiltonian <[XJ> in the form 

# = #(c) + v M +v M +v « ; (15) 

where 

H (c) =H (0) +v (c) +v ^) +v (c) (16) 

Eq. {T5J looks like an ordinary Hamiltonian with short-range potentials. Therefore the 
Faddeev procedure is readily applicable and the Faddeev components are defined by 

|V/„>=G (c '(£)v4 5) |f), a = 1,2,3, (17) 

where G^Z?) = (E — ff( c )) -1 . Since v$ is of short-range type, its acting on |f) suppresses 
those asymptotic parts in If) which are related to situations where particles jS and 7 are far 
away. With this definition of the asymptotic components, the Faddeev integral equations 
take the form 





= Gf } (E)vf' 


'(|V52} + |¥3» 


(18) 


V2) 




'(IVi}+IW) 


(19) 


¥3) 




\\Wi}+\¥2)) , 


(20) 



where 

G$(E) = (E-H® -vifV 1 . (21) 

With this treatment of the confining potential we achieved that the Faddeev decompo- 
sition acts as an asymptotic filtering, and the Faddeev integral equations formally contains 
only short-range potentials. All the long range terms are kept in the Green's operators. How- 
ever, we should be careful here. The Hamiltonian H^ c > contains the long-range confining 
parts of the original Hamiltonian. Thus H^ c ' has infinitely many discrete eigenstates. At 
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those discrete energies (E) is singular. So, acting on a vanishing v„ \W), they may pro- 
duce a non-vanishing \f/ a ) . Thus it may happen that at certain energy the Faddeev equations 
have solutions where none of the \ \ff a ) components are zero but their sum, the total wave 
function If), vanishes. These solutions are spurious or ghost solutions. They are due to the 
pole of G( C \E). However, G^(E) is just an auxiliary operator. It depends on the splitting 
of v a . Therefore the spurious solutions depends on the splitting, while the real solutions 
should not. If we vary the way we split the potential in Eq. ( |14[ l we can shift the energy of 
the spurious solutions and with a careful choice of v' c ' we can push them out of the energy 
region of physical interest. 

A further advantage of the Faddeev method that the identity of particles greatly sim- 
plifies the equations. If particles 1, 2, and 3 are identical, then in its natural Jacobi 
coordinate system {xi,yi}, looks like ifc in its natural Jacobi coordinate system {x2,y2} 
and y/2 in its natural Jacobi coordinate system {x3,y 3 }. On the other hand, the wave func- 
tion is either symmetric (p = 1) or antisymmetric (p = —1) with respect to interchanging 
two particles 

&23\Yi)=p\Vi)- (22) 
Building this information into the formalism we arrive at a single integral equation 

|V/ 1 >=2G ( 1 c) v ( / ) ^ 12 3|^i), (23) 

where £^123 = ^12^23 is the operator for cyclic permutation of all three particles ^123 1 Vi) = 
1//2). The three-quark system is completely antisymmetric in terms of color quantum num- 
bers, so it should be symmetric in the other quantum numbers. To achieve the proper symme- 
try for all the three quarks we only need to ensure a symmetric wave function with respect to 
exchange only two particles. This can easily be achieved by selecting the subsystem angular 
momentum I, spin s and isospin t such that 

(-)'+*+' = 1. (24) 

This will ensure the correct symmetry for all the three quarks in the baryon. 



4 Solution of the Faddeev equations for confining potentials 

We solve the Faddeev equations be using the Coulomb-Sturmian (CS) separable expansion 
approach. In configuration space the CS functions are defined by 



^ w/;fc ) = y {n + 21 + i) t ( 2fcr ) /+1 ex P(~M (2br) , (25) 

where n = 0, 1 . .., / is the angular momentum, b is a parameter and L is the Laguerre poly- 
nomial. The biorthonormal partner to \nl\b) is defined as (r\nl;b) = 1/r- (r\nl;b). The or- 
thogonality and completeness, in angular momentum subspace, is given by 

(n'l;b\nhb) = (nf^b\nl;b) = 8„,„ (26) 

and 

N N 




6 



The CS functions also have a nice analytic form in momentum space 

W*) = A ^ J n+ \X l)l T£ M ^ (44) (28) 
u ; yn{n + 2l + l)\ (p 2 +b 2 ) 2l + 2 " \p 2 + b 2 J 



and 

p z + b 2 
2b{n + l + l) 

The three-body Hilbert space is the direct sum of two-body Hilbert spaces associated 
with coordinates \ a and y a . Then the appropriate basis is given by 



^ = 2bin + l + l) ^- (29) 



\nlvk) a = {\nl;b x ) a ® |vA;^)„} , (30) 

where \nl;b x ) a and |vA;fc v ) a are associated with coordinates x a and y«, respectively. We 
can have three different bases of this kind for a = 1 , 2, 3, and the completeness relation takes 
the form 

N 

1 = lim Y\nlvX) a a {nlvX\ 

(31) 

= lim V |«/vA) a a (nlvX\ = lim 1^ . 

In order to have a convenient representation for the Faddeev equations, we plug the unity 
operators into the equations 

Wa) = Um G { <\e) [l N J£l»\ ¥p )+l N J°hy\ ¥Y )] . (32) 

This representation becomes an approximation if we keep N finite. This is equivalent of 
approximating the short-range potentials in the three-body Hilbert space basis 

v£ } « f \nWX) a v^ p {n'f^X'\, (33) 

n,v,n'v' 

where = a (nv\ v„' \n'v')p. This approximation turns the Faddeev equations into a matrix 
equation 

V/ a = GL e) (£)(v^^+v^ r ) (34) 

(35) 

where 

G^ ] (E) = a (nlvX\(E-H^-v^)- 1 \n'l'v'X') a . (36) 

By rearranging, we get a homogeneous linear algebraic equation for the component vector 

(&T l -iOV = 0, (37) 

where 



/(G ( i e) ) _1 (£) Q \ 

o (G^yHE) o 

V (G^VW 



(38) 
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f 



( v« vif 
$ Q vg 
VjM 



(39) 



and 



¥ 



(40) 



This equation is solvable if and only if the determinant vanishes 

det(^ 1 -r)=0. (41) 

(s) 

While the numerical evaluation of the matrix elements v a „ is straightforward by using the 



(c) 

transformation of Jacobi coordinates, to determine the matrix elements (T a ' we need further 
approximations. 



4.1 The CS matrix elements of 

(c) 

The Green's operator G a is the resolvent of a complicated three-body operator 

^ = + i + $ =H®+ vf + V f + v a . (42) 

Here 

<=L^W (43) 

is the relativistic kinetic energy operator, m, is the mass of the constituent quarks and 
k, are the individual quark three-momenta in the frame where the total three-momentum 
P = £? =1 k, = 0. Since this operator is too complicated for a straightforward evaluation of 
its resolvent we adopt the following procedure. We separate off a Hamiltonian which de- 

(c) 

scribes the asymptotically most relevant part of H a ' , and which we try to handle without 
approximations. The rest, which is not dominant asymptotically, can be approximated on 
the CS basis. 

We define the operator H a as a sum of two-body operators acting on x a and y a , respec- 
tively, 

H a =h Xa (x a )+h, a {y a ), (44) 

where 

hx a {x a ) = \jpl a +Vx a + Va (x a ) (45) 

and 

h ya (ya) = ^pl+^ a + U { a ) (y a )- (46) 

Here ji Xa = mpm r / (mp + m r ) and jj. ya = m a (mp + my) / (m a + mp + m r ) are the reduced 
masses. The operator h Xa is the semi-relativistic Hamiltonian of the (jS , y) pair, while the 
operator h ya describes the asymptotic motion of the particle a and the pair (/3 , y) . The aux- 

(c) 

iliary confining potential w a (y a ) should behave as v c „ (xp ) + v c 7 (x r ) for large y a distances. 
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(c) 

The Green's operator G a satisfies the resolvent relation 



Ga — G a +G a U a Ga , (47) 
where G a (z) = (z — H a )~ l is the resolvent of H a and 

Ua —Ha ' — Ha 

= E \/ k f + m f - \Jp 2 x a +tf a - \JPy a +rf a + vj° {xp ) + vf {Xy) - utf (y a ). ^ 
i—l 

(c) 

In calculating the CS matrix elements of G a we again make the separable approxima- 
tion in Eq. ( |47[ l and get 

(&(E)y l = (d a (E)y 1 -U a , (49) 

where G a (E) = a {nv\G a {E)\n'v') a and U_ a = a {nv\U a \n'v') a . The later matrix elements 
can also be evaluated numerically partly in configuration, partly in momentum space. 



4.2 The matrix elements G a 

For calculating the matrix elements of G a we can adopt the technique used before in the 
Coulomb case \ \2\. The method is based on the Riesz-Dunford functional calculus. Given 
a self-adjoint operator h on a Hilbert space and an analytic function / in the neighborhood 
of o{h), the spectrum of h, the function of the operator f(h) is defined by a Cauchy-type 
contour integral 

f(h) = ^.^f(z)(z-h)- 1 dz, (50) 

where c € the boundary of a domain D D o(h). 

The Green's operator G a is a function of the self-adjoint operator h Xa . Therefore 

G a (E) = (E-h ya -h Xa )- 1 

= 2 J[ j£, Sy a ( E ~ z) 8x a (z) dz , 

where g ya (z) = (z — ft ya ) _ 1 and g Xa (z) = (z — A ta ) _ 1 ■ The contour ^ should be taken such 
that it encircles the singularities of g Xa and avoids the singularities of g ya . 

In fact, both h Xa and h ya are two-body Hamiltonians with confining potentials having 
discrete spectra. Figure [T] shows the analytic structure of the integrand in Eq. \5l) . For the 
better visibility, we display the integrand for complex energy E = E r + ie, e > 0. Then the 
poles of g ya become well separated and we can easily draw a contour that encircles the 
poles of g Xa and avoids the poles of g ya . In Figure [2] the contour ^ is deformed such that 
it shrinks to a few lowest-lying poles of g Xa and opens up to an integral along imaginary 
direction. Now, even in the e — > limit, the mathematical requirements for the contour 
integral are satisfied (Figure |3j. 



Fig. 1 Analytic structure of the integrand in Eq. pi) for complex E = E r +ie. The contour V circumvents 
the poles of g Xa without incorporating the poles of g y . 



<$ <•> <•> ••• 



Fig. 2 The contour in FigureMJis modified such that it shrinks to a few low-lying poles of g Xi 
to a contour along imaginary direction. 



and opens up 



E->0 

x <•) 



X ••• X •. X i-X • X 



♦ • ♦ • 



Fig. 3 Now the contour avoids the singularities of g y even in the e — > limit. 



4.3 Matrix elements of two-body Green's operators 

To evaluate G a by a contour integral of two body Green's operators, we need the CS matrix 
elements of g Xa and g ya . Both h Xa and h Ya are semirelativistic Hamiltonians with short-range 
term 

h = sjp 2 + m 2 + air + v is \ (52) 

We can approximate the short-range potential v^ s ' in a separable way using a double CS 
basis 1131 . It amounts to expanding v( s ' as 

N . . _] . . _] 

vW w £ \nl;b) ( (nl;b\ml\b')\ {ml;b'\v^\m'l';b') Um'l';b'\n'l r ;b)) (n'l';b\ 

nmm'n' 

N 

« Y,\nl;b)v {s) {ril';b\ . (53) 

nn' 

This is an exact representation if N goes to infinity and becomes an approximation if N is 
kept finite. The CS matrix elements of the potential have to be evaluated numerically. This 
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can easily be done in configuration or momentum space, depending how the potential is 
defined. 

For the operator g(z) = (z — h) , we can immediately write down a resolvent equation 

g(z)=g^(z)+g^(z)v^g(z), (54) 
where g( c ' (z) = (z — h^)^ 1 , and 

h {c) = \/p 2 + m 2 +a\r . (55) 
With the help of the approximation ( |5"3j l we can solve equation l |54| i, and get 

g(z) = ((g {c) (z))- l -v_^)- 1 . (56) 

What remains is the evaluation of g^ (z). 

The main advantage of using CS functions in few-body Coulombic calculations comes 
from the fact that the Coulomb Green's operator can be given in analytic form. This is 
possible because in CS basis the Coulomb Hamiltonian is an infinite symmetric tridiagonal, 
or Jacobi, matrix. We found that the resolvent of a Hamiltonian having a Jacobi form can 
be given in terms of a continued fraction [14]. This result can be generalized to polynomial 
potentials |6]. A non-relativistic kinetic energy with a polynomial potential has an infinite 
symmetric band structure. If the polynomial potential is a linear potential, the band matrix is 
penta-diagonal, if the confinement is quadratic, the band matrix is septa-diagonal. However, 
band matrices can be considered as tridiagonal matrices of block-matrices. In particular, a 
penta-diagonal band matrix can be considered as a tridiagonal matrix of 2 x 2 block matrices. 

Unfortunately, the relativistic kinetic energy operator does not have a band structure on 
CS basis. But numerical studies revealed that for large n and n' values the dominant elements 
form a penta-diagonal band matrix. So, by neglecting the asymptotically not dominant el- 
ements outside the penta-diagonal band, we can use the matrix continued fraction method 

ED 

g} c} {z) = (Z-5 i)n '5^„/7„» in / +1 C„/ +1 /„/ +1;n /)" 1 , (57) 
where C is a matrix continued fraction defined recursively by 

Q+l = (Ji+l,i+l — Ji+l,i+2Q+2Ji+2,i+l) ' (58) 

where J = z — h in CS representation, refers to 2 x 2 blocks and ri = n/2. So, g c is an 
inverse of the modified J matrix. The modification is a matrix continued fraction C, and 
it only effects the bottom right 2x2 block of J. Numerical studies showed, that although 
we neglected the asymptotically non-dominant terms in J, this construction of g^ is very 
accurate 1151 . 



5 Goldstone boson exchange model for baryons 

In this work we consider a semirelativistic quark model with Goldstone boson exchange 
quark-quark interaction. The kinetic energy operator is given in the form of Eq. (43) . The 
quark-quark interaction is derived from Goldstone-boson exchange model 



_ yconf y x 



(59) 
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where confinement potential is taken in the form 

VS mf = V + Cx a . (60) 
The chiral potential is a sum of octet 

3 7 

V« (x«) = £ V % (x a )XpXyOpO r + £ V K (x a )XpXy OpO r +V v (x a )X$X*OpO r , (61) 

F=l ' ' F=4 

and singlet 

y 2 

Va(x«)= -V^(x„)CT /j a 77 (62) 

meson-exchange terms, where a and X are the quark spin and flavor matrices, V n , Vk, V-q 
and V^i represent the form factor for %, K, rj and rj' meson exchanges, respectively. In the 
model of Ref. 1161 171 this term is taken in the form 

^M = g T y i — Aj 6 - . (63) 

47T 12mpin r { x a x a ) 

The second Yukawa term is a smeared delta potential and its coupling constant Ap is as- 
sumed to have a linear dependency on meson masses 

A F =A + Kfi F , (64) 

where Aq and K are free parameters. In this model the masses of quarks and the mesons are 
fixed parameters, m u = m d = 340 MeV, m s = 500 MeV, \l % = 139 MeV, fi K = 494 MeV, 
pL-q = 547 MeV, fi^i =958 MeV, and the meson octet-quark coupling constant was adopted 
as = 0.67. The other parameters for the model were determined by fitting manually to 

the observed baryon spectra. Excellent agreement with experiment was found with Vq = 416 
MeV and C = 2.33 fmT 2 , A Q = 2.86 fin -1 , K = 0.81 and {go/gs) 2 = 1-34 Gi- 
lt is natural to incorporate y con f into v c and V x into v^. Then, to ensure that the poles 
of G^ c > do not generate spurious solutions, we add a repulsive Gaussian term to y( c ' , which 
we subtract from v(' s ) 

v (c) = yconf + aoe -(r/r f (65) 

and 

V W =V x -a oe -^ r °'> 2 . (66) 

The parameters of the auxiliary potential have been taken as «o = 3 fm - 1 and ro = 1 fm. By 
this choice of the parameter values any bound states of H( c > are avoided below w 2 GeV 
The values of oq and ro also influence the rate of convergence but not the final (converged) 
results. 
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6 Results 

In this approach we solve the three-quark problem in a discrete Hilbert space basis repre- 
sentation. To get reliable results we have to achieve convergence in terms of basis states, i.e. 
in angular momentum channels, in terms of N of Eqs. l |33| >, ( |49[ l and ( |53|l, and in terms of 
n, which denotes the highest index in matrix continued fractions. Tables[T] and [2] show the 
convergence of the nucleon and excited nucleon state, respectively, in terms of N and n. We 
used 10 angular momentum channels as specified in Table [3] We can see that N = 20 and 
n = 60 provides very accurate results. In calculating other results we kept these basis sizes 
fixed. 

Table [3] shows the convergence of the nucleon and excited nucleon states in terms of 
angular momentum channels. The quantum numbers of the channels are selected to meet 
the angular momentum algebra and the symmetry condition l |24[ >. In Table [4] we present the 
results of our detailed calculations for various light baryon states. For comparison we also 
provide the results of a variational calculation [7| and the results of our previous Faddeev 
calculations where the confinement were treated in an approximate way Q. 

7 Summary 

In this work we presented an improvement over our previous semirelativistic Faddeev cal- 
culations of light baryons [5 |. The main improvement is the exact treatment of the asymp- 
totic confinement. For this purpose we calculated the channel Green's operator by a contour 
integral of two-body Green's operators. The two-body Green's operators contain a linear 
confinement term, and they were evaluated with the help of a matrix continued fraction. 
Numerical studies show some improvement over the previous solution method, but the dif- 
ference is little, well within the range of experimental errors. So, within the framework 
of Faddeev method, we have a good numerical solution of the semirelativistic three-body 
Hamiltonians of baryons. 

Table 1 Convergence of the nucleon (N(939)) and excited nucleon (N(1710)) mass in terms of n of , the 
terms used in evaluating the matrix continued fraction. We used N = 24 is the separable expansion of the 
potential operators in the Faddeev equations. 



N 


n 


N(939) (MeV) 


N(1710) (MeV) 


24 


30 


938.08 


1775.85 


24 


40 


938.06 


1775.84 


24 


50 


938.05 


1775.84 


24 


60 


938.04 


1775.84 


24 


70 


938.04 


1775.84 


24 


80 


938.04 


1775.83 


24 


90 


938.04 


1775.83 
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Table 2 Convergence of the nucleon (N(939)) and excited nucleon (N(1710)) mass in terms of N, the terms 
used in the separable expansion of the potential operators in the Faddeev equations. The matrix continued 
fractions were evaluated up to n = 60. 



N 


n 


N(939) (MeV) 


N(1710) (MeV) 


16 


60 


938.16 


1775.86 


17 


60 


938.14 


1775.85 


18 


60 


938.12 


1775.85 


19 


60 


938.11 


1775.85 


20 


60 


938.10 


1775.85 


21 


60 


938.08 


1775.84 


22 


60 


938.07 


1775.84 


23 


60 


938.06 


1775.84 


24 


60 


938.05 


1775.84 



Table 3 Masses of the nucleon ground state and the first two positive-parity excitations for the semirelativis- 
tic GBE chiral quark model. The convergence with respect to including an increasing number of angular- 
momentum states l x , L is demonstrated; s and ( denote the subsystem spin and isospin. In all cases, N = 20 
separable terms have been included into the expansion. A comparison to the results obtained with the previous 
calculations is given. 





Channels 






Nucleon 


N*(U<W) 


iV*(1710) 


# 


lx 


h 




t 








1 














1615.6 


1991.7 


2294.2 


2 








1 


1 


954.3 


1470.8 


1911.3 


3 


1 


1 


1 





952.8 


1469.0 


1907.3 


4 


1 


1 





1 


950.6 


1465.7 


1793.7 


5 


2 


2 








946.6 


1463.1 


1792.7 


6 


2 


2 


1 


1 


939.4 


1459.0 


1782.6 


7 


3 


3 


1 





939.3 


1458.8 


1780.0 


8 


3 


3 





1 


939.2 


1458.6 


1776.9 


9 


4 


4 








938.7 


1458.3 


1776.9 


10 


4 


4 


1 


1 


938.1 


1457.8 


1775.9 


Ref. (5] 


939 


1459 


1775 




Ref. |7| 






939 


1459 


1776 



Table 4 Masses of light baryons. 



Baryon 


Ref.0 


Ref. El 


This work 


Experiment 


N(939) 


938 


939 


938.0 


938-940 


N(1440) 


1458 


1459 


1457.8 


1420-1470 


N(1520) 


1517 


1519 


1517.0 


1515-1525 


N(1675) 


1644 


1647 


1644.4 


1670-1680 


N(1710) 


1776 


1776 


1775.8 


1680-1740 


4(1232) 


1239 


1240 


1238.8 


1231-1233 


4(1600) 


1717 


1718 


1717.2 


1550-1700 


4(1620) 


1639 


1642 


1638.9 


1600-1660 
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